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We discuss entanglement entropy of gapped ground states in different dimensions, obtained on 
partitioning space into two regions. For trivial phases without topological order, we argue that 
the entanglement entropy may be obtained by integrating an 'entropy density' over the partition 
boundary that admits a gradient expansion in the curvature of the boundary. This constrains 
the expansion of entanglement entropy as a function of system size, and points to an even-odd 
dependence on dimensionality. For example, in contrast to the familiar result in two dimensions, a 
size independent constant contribution to the entanglement entropy can appear for trivial phases in 
any odd spatial dimension. We then discuss phases with topological entanglement entropy (TEE) 
that cannot be obtained by adding local contributions. We find that in three dimensions there is 
just one type of TEE, as in two dimensions, that depends linearly on the number of connected 
components of the boundary (the 'zeroth Betti number'). In D > 3 dimensions, new types of 
TEE appear which depend on the higher Betti numbers of the boundary manifold. We construct 
generalized toric code models that exhibit these TEEs and discuss ways to extract TEE in D > 3. 

PACS numbers: 



I. INTRODUCTION 

In recent years, surprising connections have emerged 
between error correction of quantum information and 
topological condensed matter phases^. At the same 
time, ideas from quantum information have proved useful 
in defining topological phases. Two dimensional phases 
with topological order, such as those realized in the con- 
text of the Fractional Quantum Hall effect, arc gapped 
phases for which the ground state degeneracy depends on 
the genus of the space on which they are defined. Recent 
work has shown that they can be identified by the entan- 
glement properties of their ground state wavefunction^r— . 
The entanglement entropy of a region with a smooth 
boundary of length L takes the form Sa = ol\L — &o7, 
where 7 is the topological entanglement entropy, bo is 
the number of connected components components of the 
boundary of region A, and we have dropped the sublead- 
ing terms. In gapped phases without topological order, 
such as band insulators, 7 = for a smooth boundary. 
These predictions have been verified in the context of 
a number of specific D = 2 models^— and in D = 3 Z2 
toric code models^. In this paper, we discuss the general 
structure of entanglement entropy for gapped topological 
and non-topological phases in D > 3. 

One notices that in D = 2, the topological entan- 
glement entropy depends only on a topological prop- 
erty of the boundary- in this case the number of con- 
nected components. There are two equivalent ways of ex- 
tracting the topological entanglement entropy^. First, 
via the scaling of the entropy with boundary size for 
smooth boundaries, to extract the constant term. The 
second, by considering a combination of entanglement 
entropies of three suitably chosen regions A, B, C, so 
that —7 = Sa + Sb + Sc — Sab — Sbc — Sac + Sabc- 

Here, we will discuss analogous questions in D > 3. 
In particular, consider a gapped D — 3 phase, and a re- 
gion A with a smooth boundary. (1) If the entanglement 



entropy Sa contains a constant term, does it necessarily 
reflect topological order? (2) The boundary of A is a 
closed two dimensional surface that has two topological 
invariants associated with it - the number of connected 
components, and the genus (number of handles) of each 
component. Does this imply there are two distinct types 
of TEEs and correspondingly two varieties of topological 
order in D = 3? The answer is no to both these ques- 
tions, as we elaborate in this paper. We show that even 
a trivial gapped phase, with no topological order, one 
can have a constant term in the entanglement entropy 
in D = 3 (and in any other odd dimension). Hence, 
this by itself does not signify topological order. More- 
over, this constant is generally genus dependent, ruling 
out a topological origin for a genus dependent entangle- 
ment entropy. This reduces the number of possible TEE 
to the same as D = 2. We discuss generalization of the 
Kitaev-Preskill scheme^ to extract the TEE in D = 3; 
and why some naive extrapolations fail. 

A deeper understanding of TEE is obtained by consid- 
ering higher dimensions. We show that at least one new 
topological entanglement entropy appears on going up 
every two dimensions. Thus, while D — 2 and 3 are simi- 
lar, a new topological constant docs appear in D = 4 (in 
D = 2nand2n+l, there are thus n constants). These are 
related to the Betti numbers^ of the boundary. We con- 
struct topological phases that manifest these new TEEs, 
and explicitly calculate their value. These are based on 
a discrete gauge group G. In all dimensions, the TEE 
for discrete gauge theories is — log \ G\ per connected sur- 
face component, where \G\ is the number of elements in 
G. These theories capture both abelian (like the Z2 toric 
code) and non-abelian phases and the ground state of 
these theories correspond to condensate of closed loops. 
One can also consider more general abelian discrete gauge 
theories where the fluctuating loops are readily general- 
ized to fluctuating p-dimcnsional surfaces. These man- 
ifest explicitly in the topological entanglement entropy, 



through the appearance of new topological constants that 
depend on higher Betti numbers. Furthermore, a previ- 
ously discussed duality between p and D — p theories in 
D dimensions 1 -^—, is reflected in the structure of TEE. 

To isolate topological contributions it is useful to know 
the structure of entanglement entropy in trivial gapped 
phases. Since correlations are local in such phases, we 
propose an expansion of entanglement entropy Sa based 
on adding individual contributions from patches on the 
surface of region A: Sa — Y^,i The entropy densities 
Si will depend on local properties, such as the local cur- 
vature of the surface. One can then expand the entropy 
density in polynomials of curvature and its derivatives, 
similar to the Landau expansion of free energy density 1 ^. 
The contribution from higher order terms to Sa are sub- 
dominant for large surfaces. Interestingly, not every term 
is allowed in this expansion. When we divide space into 
a region A (inside) and A (outside), the entanglement 
entropy of both are equal i.e. Sa = Sa- This imposes a 
Z2 symmetry on the expansion that is unique to ground 
state entanglement entropy^. This has important con- 
sequences. Consider for example, D = 2. The entropy 
density is not allowed to depend linearly on the bound- 
ary curvature k, which changes sign on interchanging in- 
side and outside. Thus, the expansion of entropy density 
for a trivial D = 2 phase is Si = oq + 0,2 k 2 (7^) + . . . , 
which when integrated around the boundary leads to 
Sa = ot\L + ots/L + . . . , where L is the length of the 
boundary of region A. The first term is the area law, 
and the next term is two orders of L down, due to the 
Z2 symmetry, which eliminates the constant term in to- 
tal entropy for smooth boundaries. Thus the existence 
of a constant term in a gapped D = 2 state implies a 
non-trivial phase i.e. topological order. In general, this 
method predicts that for an isotropic, parity invariant 
state without topological order, the entanglement en- 
tropy in even (odd) spatial dimensions depends only on 
odd (even) powers of L, the linear scale of the boundary. 

We will assume it is possible to take the continuum 
limit for all the phases that we consider in this paper. 
This assumption exclude phases such layered Zi 

topological ordered phases in D = 3 whose topological 
entropy depends on the local geometry of the region A. 

The paper is organized as follows: in section [TT] we dis- 
cuss the general structure of entanglement entropy for 
gapped phases and explain the basic assumptions under- 
lying our discourse. In section UTT1 we introduce the afore- 
mentioned curvature expansion for entanglement entropy 
of trivial gapped phases and study its consequences. In 
section HVl and IVl we study topological ordered phases in 
D = 3 and D > 3 respectively, through extracting the 
dependence of entanglement entropy of a region on the 
topology of its boundary. We also generalize the con- 
structions for extracting topological entrop y 4 ' 5 ' 10 . 
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Figure 1: The local part of the entropy of region A is the sum 
of contributions of small patches on the boundary. 

II. STRUCTURE OF ENTANGLEMENT 
ENTROPY FOR GAPPED PHASES 

In this article, we will assume that the entanglement 
entropy of a region A can be decomposed into two parts: 

Sa = Sa,iochi + Sa .topological (l) 

We postpone the underpinnings of this assumption to 
Sec II VI when we study topologically ordered phases. Here 
Siocai is defined by the property that it can be written as 
a sum over contributions from patches located along the 
boundary of region A: 

SA,local = ^2 Si (2) 

i 

where Si depends only on the shape of the patch i, 
and not on the rest of the surface or how it fits with 
other patches, see FigfTJ at least if the edge of the patch 
connects smoothly to all other patches. 

We assume the other contribution St opo iogicai is topo- 
logically invariant, i.e., it does not change as the bound- 
ary is deformed unless the topology of the region changes. 
If such a term is present and if it cannot be expressed in 
a local way, then the phase has long range entanglement, 
which is the hallmark of topological order^ 4 -! 5 -. 

Let us consider the assumptions under which the de- 
composition (Eqn(2]) would be possible for a trivial (i.e. 
not topologically ordered) gapped phase. The reduced 
density matrix corresponding to a region A for the ground 
state wavefunction may be written as 

Pa = e~ HA /Z 

where Ha is the so called entanglement Hamiltonian and 
Z = tr(e~ HA ) so that tr(pA) = L Therefore, we can 
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think of pa as the thermal density matrix at temperature 
T = 1 for the Hamiltonian Ha and the von Neumann 
entropy Sa = —tr(pA logpA) as the thermal entropy for 
this system. Let us define Pa(T) = e~ HA l T /Z(T) where 
Z(T) = tr{e- HA ' T ). Clearly, p A = P~a(T = 1) and S A 
obeys the following equation: 

9=9 - rdTd(H A ) 

— dA,local — J -j, q^, \d) 

where (Ha) denotes the thermal average of Ha at tem- 
perature T with respect to the density matrix p(T). We 
claim that the entanglement entropy would admit an ex- 
pansion such as EqnJ2]if the following conditions are sat- 
isfied. 

• Ha can be written as a sum of local operators O's 
i.e. H A = E x O(x). 

• There is no phase transition for the Hamiltonian 
H A for T > 1. 

The first condition along with the fact that all corre- 
lations of local operators are short-ranged in a gapped 
phase imply that Ha has non-zero support only near the 
boundary of region A (within the distance of correlation 
length). In other words, the degrees of freedom inside 
and outside of region A are coupled only through oper- 
ators that lie within a distance ~ £ from the boundary. 
This implies that (Ha) — X)i ^« where i denotes a point 
at the boundary of region A and the hi's depend solely on 
the properties of the boundary in the vicinity of point i. 
The second condition implies that the integral in EqnJ3] 
does not admit any singularity so that all terms Si in the 
EqnJ5] arc finite. Physically, this means that the actual 
system of interest is smoothly connected to its T = oo 
zero correlation length system where EqnJ5] holds triv- 
ially. 

III. ENTANGLEMENT ENTROPY OF TRIVIAL 
GAPPED PHASES 

In this section we will focus on understanding the 
leading and sublcading dependence of Sa,Ioco.i on L for 
gapped trivial phases of matter. Let us assume that the 
boundary of region A is smooth, and further that the 
phase is isotropic and parity invariant (consequences of 
the violation of these assumptions are discussed at the 
end of this section and in Appendix [B]). Then, as we 
will see below, in the absence of topological order, only 
alternate terms in the power series expansion of Sa(L) 
appear: 

S A ,local(L) = ai L°- 1 + a 3 L°- 3 + a 5 L D - 5 + ... (4) 

i.e. only those with odd co-dimension exponent can ap- 
pear. This expansion implies a distinction between even 




S(A) = S(B) 

Figure 2: Illustration of the Z2 symmetry for the curvature 
expansion discussed in the text. 

and odd dimensions: in even dimensions, any constant 
contribution to the entanglement entropy must come 
from SA,topoiogicai, and thus indicates topological order. 
In odd spatial dimensions, a constant term may appear 
in the local entropy, making it more difficult, though still 
possible, to isolate topological contributions (note that 
these conclusions apply only to smooth boundaries in ro- 
tationally symmetric systems; a corner can produce a 
constant term for even dimensions as well as odd ones). 

Let us consider some instances of Eqn0]that will moti- 
vate its derivation. First, it is well known that in D = 2, 
a constant term in Sa(L) implies the presence of topo- 
logical order— On the other hand, in D = 3, a constant 
term can appear for a non-topologically ordered phase, 
such as a gapped scalar field. Consider for a moment^ a 
massless field, where it is known that the entanglement 
scales as Sa ~ L 2 + log(L) for a spherical ball of radius 
L. Now, providing a mass m to the scalar field would 
cut off the log(L) term and instead lead to a constant 
contribution proportional to log(l/m). We verified this 
explicitly using a numerical calculation similar to RefJ^. 
Interestingly, when the surface of region A is flat, then 
there is no constant contribution as we show in Appendix 
A (a flat boundary is possible when the total system has 
the topology of a three torus T 3 ; then A may be taken 
to have the geometry T 2 x I, where I is a line segment). 
This indicates that the presence or absence of a constant 
term may have something to do with the curvature of the 
boundary of region A. In the next subsection, we make 
this statement precise and explain the observations made 
above. 



A. Entropy Density Functional 

Let us consider a region in two dimensions for con- 
creteness. We postulate that the local entropy Sa.Ioco,i is 
given by the following integral: 

S A = Y,Si= j daF(K,d K ,...) (5) 
where F(n, <9k, . . . ) is the "entropy density functional" . 



4 





Full symmetry 


All symmetries broken 


Broken parity alone 


D = 2 


All j T () T —~'Z 

All even terms L , L 


Constant term 


All even terms 


D > 3 odd 


All Odd terms L u ~'\ L"" 4 , . . . 


Nothing forbidden 


Odd terms with positive ex- 
ponents and also L _1 if D = 
3(mod 4) 


D > 4 even 


Even terms L"-' 1 , L"- 1 , . . . 


Nothing forbidden 


Even terms 



Table I: Terms in the entropy forbidden by symmetries. The three columns describe systems with rotational and parity 
symmetry, no spatial symmetry, and rotational but not parity symmetry. The entries list the scaling of terms L k that are 
forbidden from appearing in the entropy. 



In a gapped phase, the entropy Si of a patch larger than 
the correlation length can depend only on the properties 
of the patch, such as its length Aer^ and curvature k^, as 
well as derivatives of the latter, and must be proportional 
to A(Ti. Hence Si = AaiF(Ki,d n Ki). Taking the limit 
where the patches become microscopic compared to L 
(but greater than £) leads to Eq. ([5]). 

The entropy density functional always satisfies a Z 2 
symmetry, which is the key to understanding the In- 
dependence of the entropy. The symmetry results from 
the fact that, if A and B are complementary regions, 
then Sa = Sb- Therefore, changing 'inside' to 'outside' 
keeps the entanglement entropy invariant. Now, under 
this transformation, radii of curvature clearly change sign 
k — > —k, and this constrains the entropy density func- 
tional (Fig(2|). As an illustration of this Z 2 symmetry, 
consider the form of the functional F for a gapped two 
dimensional system. On a smooth boundary, one can ex- 
pand the function F in a Taylor series, retaining the first 
few terms: 

<9™k;) = a + aiKi + a 2 n 2 + b 2 d a Ki + . . . (6) 

The first term gives the boundary law Sa = a^L. The 
second term, if it is present, would give a constant con- 
tribution § dan = 2tt for the curve shown, which would 
be a non-universal constant contribution. However, such 
a term is in fact forbidden by the Z 2 symmetry, since 
the term is odd in k. The term k 2 gives the next con- 
tribution to the entropy that is proportional to § dan 2 . 
If the shape of region A is kept fixed, then this contri- 
bution scales as 1/L. The term d a n is allowed by the 
Z 2 symmetry (since both the derivative and the radius 
of curvature change sign, assuming that the direction of 
the curve is set by a 'right hand rule', whereby the arc 
length increases along a specific direction) , yet it still van- 
ishes because it is a total derivative. Generalizing these 
arguments, one finds that Sa — YlkLo a ik+iL l ~ 2k . 

In general dimensions, we find similar results if we con- 
tinue to assume rotational, parity, and translational sym- 
metry. These assumptions imply that F can depend only 
on the metric tensor g a/3 and on the extrinsic curvature of 
the surface. The latter tensor does not appear when con- 
sidering intrinsic properties of a manifold (as in general 
relativity). However, entanglement entropy does depend 
on the embedding of the boundary dA, since it is mea- 
suring the entanglement of the degrees of freedom in the 



space around the surface. The extrinsic curvature is a 
tensor k Q( 3 with two indices (see Appendix [F] for a short 
primer on the requisite differential geometry). Thus each 
term in F contains some number of factors of K a p and 
its covariant derivatives, with all the indices contracted 
by factors of g lS (if parity is broken, the antisymmetric 
volume tensor -y Q i - Q «-i i s allowed as well, Appendix IB1. 

The inside/outside symmetry further limits the form of 
the terms in F: it implies that each term in F includes 
an even number of factors n K of k. The total order of all 
the derivatives nr> must also be even. This follows from 
rotational symmetry. For rotational symmetry to be re- 
spected, one has to contract all the lower indices with 
the tensor g lS . This leads to an even number of lower 
indices, that include the derivatives as well as the curva- 
ture indices. Since the curvature tensor is of even rank, 
the number of derivatives no has to be even. Putting ev- 
erything together, one finds that the contribution to the 
entropy density F scales as that clearly has 

an even exponent, explaining why only alternate terms 
appear in the entropy, Eq. [U 

When rotational or inversion symmetries are broken 
spontaneously or by applying a field, additional terms 
appear in the entropy, as summarized in Table [I] We 
provide the details leading to these results in Appendix 

m 

The local entropy can also contain topology-dependent 
terms e.g. the term J GdA = Attx in three dimensions 
where G is the Gaussian curvature which is the determi- 
nant of the matrix K a p . Note that this term is compatible 
with the symmetry k — > —k, since G is quadratic in k. 
Hence, as mentioned earlier, the presence of a term in the 
entropy that is proportional to the Eulcr characteristic 
does not necessarily correspond to topological order. In 
general dimensions, J detKGL4 is topological, but it is only 
symmetric in odd dimensions, where it is proportional to 
the Euler characteristic of the boundary (in general, it 
is proportional to the Eulcr characteristic of the region 
itself^). 



IV. TOPOLOGICAL ENTANGLEMENT 
ENTROPY IN D = 3 

We now turn to the topological part of the entangle- 
ment entropy. Our starting point is Eqn. [T] which we 
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rewrite here for convenience: 

local + S A ^topological (7) 

This decomposition is what enables the extraction of 
topological entropy using Kitaev-Preskill^ or Levin- Wen^ 
constructions for two-dimensional topological ordered 
phases. The assumptions underlying this equation are 
somewhat tricky. Though Eqn(7] holds for toric code 
models in all dimensions and there is strong numerical 
evidence that it also holds for many interesting two di- 
mensional topological ordered states such as Z2 spin liq- 
uids, quantum dimer models and various quantum Hall 
stateaSASi, Eqn.7 surreptitiously rules out a layered Z2 
topologically ordered state. Such a state would lead to 
a correction in entanglement entropy ASa = —"f2DL z , 
for layering perpendicular to the z direction. Here 72D 
is the topological entanglement entropy associated with 
the theory living in each layer. Clearly, AS a is not topo- 
logically invariant. The assumptions underlying Eqn. 
most likely also do not apply to the self-correcting code 
state of Refi22; this state (whose ground state degeneracy 
depends on the divisibility of system size by powers of 2, 
for example) illustrates why we need to make an assump- 
tion of this type. Nevertheless, we will briefly discuss the 
topological entanglement entropy of layered Z2 state in 
Appendix [DJ 

Independent contributions to St op ological '■ The bound- 
ary dA of a three dimensional region A is a compact 
manifold that is characterized by Betti numbers 60 and 
b\ (note that for compact manifolds 62 = bo). As we 
show in Appendix IE1 in three dimensions S ' A,topoiogicai is 
a linear function of bo, &i, say, S t0 p iogicai = -7060 - 71&1 
(we assume that the space in which region A is embedded 
has the topology of R 3 , otherwise more complicated de- 
pendence is possible in principle) . This might lead one to 
suspect that there are two different kinds of topological 
orders in three dimensions, namely, those corresponding 
to a non-zero 70 and 71 respectively. However, 60,61 are 
related to the Euler characteristic x through 26o — 61 = X- 
Thus, one may redefine S' Alocal = Sa,Ioco,1 + "X an d 

^^topological = & A 

.topological ®QC 

without changing the 
entropy and a may be adjusted so that the b\ depen- 
dence of Stopological is canceled out. Here the term ax 
may be thought of as both local and topological. It is 
local, because the Euler's formula, x = V — E + F gives 
a local expression for this term, where V, E, and F are 
the number of vertices, edges, and faces into which dA 
is divided (alternatively, in a continuum theory, ax can 
be incorporated into the entropy density F since x is 
the integral of the Gaussian curvature). It is also topo- 
logical, because ax is independent of how the surface is 
divided up into regions. The upshot of this discussion is 
that there is only one kind of topological entropy in three 
dimensions. 

Zi string and Z2 membrane models: As an alternative 
way to understand the above result, let us study specific 
models whose topological entanglement potentially de- 
pends on different Betti numbers. Consider the following 



model of Z2 gauge theories consisting of spin- 1/2 degrees 
of freedom that live on the links of a three-dimensional 
cubic lattice: 

H r string = ~ ^2 J] T *,l ~ h T *,l ( 8 ) 

□ lea i 

where □ denotes a plaquette of the cubic lattice, the 
operators t x i,t x i live on the links I of the lattice. The 
above Hamiltonian is supplemented with the constraint 
('Gauss law') ]1/ G vertex T xJ = 1 to impose the absence of 
Z2 charges in the theory. Because of this constraint the 
gauge invariant degrees of freedom in this model consist 
of closed loops C on the edges of the lattice. In the de- 
confined phase of the gauge theory, \h\ <C 1, the loops 
condense because they do not cost much energy. The en- 
tanglement entropy of this model for a region A depends 
only on the Betti number 60 of dA, since each compo- 
nent of the boundary places a separate constraint on the 
loops that intersect the boundary dA. Let us take Ki- 
taev's 'toric code limit' of the above model^ by setting 
h = 0. In this limit, the constraint commutes with the 
Hamiltonian and can be included as a part of it. Hence 
the model may be written as 

H s tring.h=0 = ~ ^ IT 7 " 2 '' - Tx ' 1 
□ /£□ vertex vertexGl 

(9) 

Interestingly, the ground state of EqnJ5] may be rein- 
terpreted as a superposition of closed membranes. This 
is seen as follows. The first term in the Hamiltonian may 
be regarded as the constraint ILeD 1 " 2 ^ ~ 1- Now con- 
sider the dual lattice each of whose plaquettes is pierced 
by a link T of the original cubic lattice. A surface can 
be defined by the plaquettes of the dual lattice pierced 
by t Zj i = — 1 bonds. Due to the constraint, this surface is 
closed. Thus there is no distinction between condensed 
loops and condensed membranes in this cas o 13 ' 15 , consis- 
tent with the fact that in three dimensions there is only 
one kind of topological entanglement entropy. 

Discrete gauge theories in D = 3: Before moving on 
to the discussion of topological entropy in general dimen- 
sions, let us derive the entanglement entropy correspond- 
ing to a discrete gauge theory with general gauge group 
G for Kitaev model^ on a cubic lattice^ 3 - 

H = -t^2 %i<?2.9354 = e)-Vj2 L l L2 g L g L 4 g (10) 

Here 'p' stands for a plaquette, 's' for a star (i.e. six 
links emanating from a vertex) while <?'s are the elements 
of group G with size of group being \G\. For non-abelian 
groups one needs to chose an orientation of the links so 
that for opposite orientations, the group element on a 
link is g and The operators L a live on the links 
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and and their action is described by L gi \g 2 ) = \g1g2) or 
L gi \g2) — 1 52<7iT X ) depending on whether gi points away 
from or towards the vertex at which the action of L g is 
being considered. The ground state of |<&) of H is given 
by 



i*>= E 

{9} j 9i 92 939 a — e V plaquettes 



(ii) 



Let us divide the entire system into region A and B 
and assume that the boundary is made up of plaquettes 
of the lattice. The links along the boundary are labeled 
by the group elements h\,li2, ■•■,h n - The Schmidt de- 
composition of |$) reads 



|*> = E 

{h} 



Ah] 



Ah} 
f out 



(12) 



where 



Ah} 



E 



\{g}) 



(13) 



— e V plaquettes^: A, 
gi—hi for i£dA 



and \<f>)iut is defined similarly. All the states in the 
Schmidt decomposition enter with the same weight and 
are orthogonal, therefore the entanglement entropy is 
the logarithm of the number of states. These may be 
counted by finding all the configurations for {h} that 
satisfy the following constraint: the product of the {h}'s 
around any closed loop on the boundary must equal the 
identity's. This includes contractible as well as noncon- 
tractible loops on the surface, and each independent loop 
reduces the total number of configurations by a factor of 
\G\, leading to 



S = log{\G\ v - 1 ) 
= Vlog(\G\)- 
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(14) 



where V is the number of vertices on the boundary 
and 7 = log{\G\) is the topological entanglement entropy. 
This result for topological entanglement entropy is iden- 
tical to that for discrete gauge theories in D = 2. 





(a) 



A 






C 


B 






(b) 




(c) 



Figure 3: Fig. (a) and (b) show two valid ABC constructions 
(Eq. I15|l in three dimensions that can be used to extract the 
topological entanglement entropy. In Fig. (a) the cross-section 
of a torus has been divided into three tori A, B and C while 
in Fig.(b) a torus that has been divided into three cylinders 
A,B and C. The Fig.(c) shows an invalid construction as 
explained in the text. In all three figures, we define region D 
to be the rest of the system. 



extension of the Levin- Wen scheme to D = 3 ) . The con- 
struction involves three regions A, B, C embedded inside 
region D: 



—Itopo = Sa + Sb+Sc — Sab — Sbc — Sca+Sabc (15) 

In two dimensions, the regions are taken to be three 
120° segments of a circle. In three dimensions, a direct 
generalization of the two-dimensional construction (di- 
viding a cylinder into three sectors as in FigJSt) fails to 
be topologically invariant because the changes in the en- 
tropy near the points at the top and bottom of the cylin- 
der where A, B, C and D all meet do not cancel. How- 
ever, the regions such as the two shown in FigOk and 3b 
can be used where such points do not exist. For exam- 
ple, if one deforms the circle at which regions A, B, D all 
meet (in cither geometry), then 



A. Extracting Topological Entanglement Entropy 

in D = 3 

In the spirit of Ref4^, we would like to combine the 
total entanglement entropies of certain regions in such a 
way that the local part of the entropy cancels out while 
topological part survives. 

The Kitaev-Preskill construction, which succeeds in 
this task in two dimensions can be modified so that it 
works for three dimensions as well (RefJ^ describes an 



A7 topo = -[A(5a 
= 



S CA ) + A(S B - S BC ) + A(S A bc 



The last equation follows because each of the three 
terms in the brackets could be thought of as the dif- 
ference between entropies of two regions that differ by 
addition of region C (that is located far from the point 
where A, B and D meet). Since each region has a single 
boundary component, "f top o = S 'topological- In Appendix 



Sab)] 
(16) 
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[Cl we detail the general requirements for a construction 
that would always yield a topological invariant. 

Based on our earlier discussion of curvature expansion 
for entanglement entropy, we note that in the special case 
of a completely flat boundary between a region and the 
rest of the system, the constant term in the entangle- 
ment entropy corresponding to that region can indeed be 
identified with topological entanglement entropy 2 ^. This 
can be realized by taking the total system to be T 3 and 
region A as T 2 x I where I is a line segment (similar to 
the calculation of entanglement entropy for a free scalar 
in Appendix [AJ . 



V. TOPOLOGICAL ENTANGLEMENT 
ENTROPY IN D > 3 

Independent terms in Stopoiogicai in arbitrary dimen- 
sions: Following our discussion of topological entangle- 
ment entropy in D = 3, in this section we study the 
independent contributions to Stopoiogicai in a general di- 
mension D > 3. The boundary dA of a _D-dimensional 
region A is a compact manifold that is characterized by 
Betti numbers, bo, . ■ . , bo-i that describe various orders 
of connectivity of the surface (see e.gJ^). 

We will assume a linear relationship, Sa = 
— X^fcLo 1 Ikbk- In principle, in higher dimensions the en- 
tanglement entropy could depend on more subtle topo- 
logical properties of the boundary, but we will focus only 
on this form. Further, as we will see below, this form 
turns out to be sufficient for Kitaev models that describe 
discrete p-form gauge theories {jp > 1) in arbitrary di- 
mensions. 

To see how many types of topological entropy can exist 
in higher dimensions, first note that for compact mani- 
folds, the Betti numbers have a symmetry, = bo-i-k 
and hence the sum may be cut short, at k = [ " P ~ 1 j . Fur- 
thermore, owing to the relation x = ^2k=oi~^-) k ^ k ' in all 
odd dimensions a part of the topological entropy may be 
absorbed into the local entropy, reducing the number of 
coefficients by one more. Hence there are n topologically 
nontrivial contributions to the entanglement entropy in 
2n and 2n + 1 dimensions: 



discrete gauge theories to understand these results. For 
example, on a hypercubic lattice in D = 4, the string 
and membrane theories describe very different ground 
states^ and unlike D = 3, the membrane theory is now 
dual to itself, not to the string phase. Explicitly, in the 
'toric code limit— ^ these two theories are given by 



Hstring = _ E II Tz ' E Tx ' 1 ^ 18 ) 

□ /£□ vertices vertexGl 



membrane 



e n^.°- e ( i9 ) 



led cubes □£ cube 



As we show now, the entanglement entropy of the 
model in Eqn. [TH]in four dimensions depends on the Betti 
number bo of dA while that corresponding to model in 
Eqn. [T9l depends on the difference b\ — bo- For the sake 
of generality, let us derive the entanglement entropy of 
a generalized toric model in arbitrary spatial dimensions 
D whose ground state is given by sum over all closed d g 
dimensional membranes. This ground state describes de- 
confined phase of a d g -form abelian gauge theory. These 
membranes intersect the boundary dA of region A in 
closed membranes of dimension d g — 1, with the restric- 
tion that these intersections are always boundaries of a 
membrane of dimension d g contained in dA. For exam- 
ple, consider the entanglement of membrane model in 
Eqnll9l in D = 3 when the boundary of region A is a 
torus T 2 (note that the form of Hamiltonian for mem- 
brane theory is identical in D = 3 and D = 4). When a 
closed membrane intersects dA = T 2 , one sees that one 
can only obtain an even number of closed loops along any 
non-contractible cycle of T 2 , which would therefore form 
the boundary of two dimensional membrane. Returning 
to the general case, let us denote the number of inde- 
pendent n-dimensional membranes that belong to dA by 
C n and those that are boundary of a n + 1-dimensional 
membrane by B n . 

Using the definition of Betti numbers^ 2 - and simple lin- 
ear algebra, one finds that the entanglement entropy Sa 



d a -l 



d a -l 



Sa 



topological 



—70^0 — 7161 — • • • — 7d_ 1 6d_ 1 , if D is 
—7060 — 71&1 + • ■ ■ — 70-3 bn-3 , if D is 

V 2 2 

(17) 

Precisely such a hierarchy of states associated with dif- 
ferent Betti numbers has been arrived at by Ref.— by 
constructing a sequence of Kitaev 'toric-code' type mod- 
els where the ground state is a superposition of all p- 
dimensional manifolds on a lattice (for 1 < p < D — 1). 
This state is dual to the superposition of all q = D — p 
dimensional manifolds, so the number of distinct models 
is L^J , the same as the number of types of topological 
entropies. 

Stopoiogicai for gauge theories in arbitrary dimensions: 
Similar to three dimensions, one may study models of 



Sa«^ R d '- 1+ "C n - {-Y h - 1+n K (20) 

n=0 n=0 

° VC §ince the C n are expressed in terms of local quantities 
°sAibh as the number of edges, vertices etc. that lie on the 
boundary without any additional constraint, we identify 
the first sum as Si oca i and the second as Stopoiogicai ■ The 
proportionality constant depends on the gauge group and 
akin to three dimensions equals log(\G\) where \G\ is the 
number of elements in the abelian gauge group (note that 
the calculation of TEE in D = 3 (Eqn. [T3]) applies to 
abelian as well as non-abelian discrete gauge theories). 
Therefore 



Stopoiogicai — —log(\G\) 2, {"Y 



-1+nj 



n=0 



(21) 
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Extracting topological entropy in Four Dimensions We 
will restrict our discussion of extracting S top oiogicai to 
four dimensions. For a given region A, from Eqn. [17] 
one has Sa = Sa,IocoX — ^o7o — &i7i and one would 
like to have a construction similar to Levin- Weni and/or 
Kitaev-Preskill^ that enables one to extract the topolog- 
ical numbers 70 and 71. 

We extract 71 by a generalization of the construction 
in^. Let region A have the topology of B 2 x S 2 . Region 
B is A with a channel cut in it and has topology of B 4 . 
Finally, region C has a second identical channel cut out 
opposite to the first one and has topology S 1 x B 3 

Now, Sa — 25b + Sc is topologically invariant just as in 
two dimensions and the Betti numbers of the bounding 
surfaces are: 



b (dA) = 1, h(dA) = l 
b {dB) = 1, bi{dB) = 
b (dC) = 1, h(dC) = l 

Hence (S A - S B ) - (S B - S c ) = Since 71 ^ for 

the membrane Kitaev model H membrane (Eq fT^)) while it 
is zero for the string model H s tring (Eq fl8| in d = 4, this 
construction measures membrane correlations. 

To isolate 70, the analogous procedure, but with A 
being B 3 x S 1 suffices. The combination (Sa — Sb) — 
(Sb — Sc) gives 70 + 71, and this may be combined with 
the previous construction to extract both 70 and 71 . This 
construction selectively measures string correlations since 
70 + 71 is zero for H membrane . 

The Kitacv-Prcskill construction of dividing a disc into 
three triangles that meet at the center is readily extended 
to any even dimension. In D = 4 consider dividing the 
ball B A into five 'pentahedra' that meet at the center. 
The combination 

A({5}) = ^ - Sij + ... + S12345 (22) 

i i<j 

is topologically invariant and gives —70. Here 5i 1 i 2 ...i„ 
denotes the entanglement entropy corresponding to the 
region A il U A i2 ... U A in . 

VI. DISCUSSION AND CONCLUSION 

In this paper, we discussed the qualitative structure 
of the entanglement entropy for gapped phases. We 
introduced the concept of 'entanglement entropy den- 
sity' whose integral over the boundary of a region A 
yields the entanglement entropy of region A. For gapped 
trivial phases the symmetry constraints on the entropy 
density, including the inside-outside exchange symmetry 
Sa = S^, naturally lead to the leading and subleading 
dependence of the entanglement entropy on the linear 
size of a given region. 



In the second half of the paper, we studied the topo- 
logical entanglement entropy Stopoiogicai of topologically 
ordered systems in various dimensions. A key result was 
that in D = 3 there is a single category of TEE, as in 
D = 2, that depends linearly on the number of connected 
components of the boundary. This constrains the possi- 
ble forms of topological order in D = 3. 

We briefly discussed TEE in higher dimensions - using 
generalized Kitaev toric code like models (i.e. deconfined 
phases of p-form discrete gauge theories) to realize vari- 
ous topologically ordered phases. In D = 4 we find two 
categories of TEE . In general, one new category of TEE 
appears each time the dimension is raised by two. This 
even-odd effect is understood as follows. Stopoiogicai de- 
pends on the Betti numbers of the boundary of region A. 
In odd spatial dimensions, the Gauss-Bonnet theorem re- 
lates Betti numbers to the curvature of the boundary of 
region A. This implies that there is one linear combina- 
tion of Betti numbers that can be expressed as an integral 
of a local property of the boundary (such as curvature), 
and is thus not an independent topological contribution 
to the entanglement entropy. We also mentioned how 
to extract Stopoiogicai by a generalization of the D = 2 
Kitaev-Prcskill and Levin- Wen constructions. 

Potentially, in D > 4, Stopoiogicai may depend not 
only on Betti numbers of the boundary manifold, but 
on more subtle topological properties such as its homo- 
topy group. If such phases do exist, then entanglement 
entropy could shed light on the classification of mani- 
folds. Lattice 3D models realize a richer variety of topo- 
logical phases than the isotropic phases considered here. 
For example, there exist layered Z2 topologically ordered 
phases, which retain a two dimensional character despite 
coupling between layers. Another example is the self cor- 
recting quantum memory of Ref^. For these, the sepa- 
ration between the topological and the local part of the 
entanglement entropy is not obvious. General statements 
about entanglement in such topological phases remain for 
future work. 

One might also consider a curvature expansion for the 
fluctuations of a conserved quantity such as particle num- 
ber or total spin, inside a region A. Intuitively, these 
would be a property of the boundary of region 
Indeed, akin to entanglement entropy, one has Fa = Fj 

where F A = ^((EreA °r) 2 ) - (T, r eA°r) 2 is th e vari- 
ance of O inside the region A. Therefore, a curvature 
expansion for Fa would inherit many of the arguments 
we used to derive the leading and sub-leading behavior 
of the quantity Fa , and can provide a framework to un- 
derstand known results^—. 

Finally, it may be possible to learn more about the sys- 
tcmatics of the size dependence of entanglement entropy 
in gapless phases by a generalization of the curvature ex- 
pansion under certain conditions. Many gapless systems 
such as massless scalar/Dirac fermion also follow an area 
law and have an expression for entropy with interesting 
parallels to EqnQJ 
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Appendix A: Absence of constant term for massive 
scalar in the absence of curvature 



We are interested in the entanglement entropy of a 
massive scalar field in three dimensions when the region 
A has a geometry T 2 x I where the torus T 2 extends along 
the directions 1,2 while / is a line segment of length / 
along the direction 3. In particular, we want to show 
that the constant part of the entanglement entropy is 
in fact exactly zero i.e. S = AL 2 + 0(1/ L). This is 
consistent with the curvature expansion of entanglement 
entropy (Eqn. O since now the region A does not have 
any extrinsic curvature. 

The Euclidean action S is given by 



2tt7 2 L 2 



-2-kYL 



tsin(2Trt) 



i/L Jt=o m2 + 7 2 (2 — cos(2nt) — cos(2iru)) 



tsin(2irt) 



2tt7 2 £ 2 



m 2 + 7 2 (l-cos(27rt)) 

tsin(2nt) 



u =o Jt=0 m2 + 7 2 (2 - cos(2irt) - cos(2iru)) 



Clearly neither I\ nor I 2 contribute to a constant term 
in the entanglement entropy. Therefore, as anticipated 
from the curvature expansion, the entanglement entropy 
S = Ii + I2 does not contain a constant term for the 
T 2 x I geometry and is proportional to L 2 upto O(L ). 



Appendix B: Additional terms in Entanglement 
entropy in the absence of rotational/parity 
symmetry 



s 



\4>(ki,k 2 ,k 3 ,uj)\ 2 (m 2 +u} 2 +^ 2 (3-cos(ki)~cos(k 2 )-cos(k 3 ))) 



Broken Rotational Symmetry 



(Al) 

where we impose periodic boundary conditions in all 
directions and we have set the lattice spacing to unity. 
Periodic boundary conditions imply momenta k\ , k 2 re- 
main good quantum numbers even after making the par- 
tition. The total entanglement entropy may therefore be 
written as 



S=J2 S 1D (M(k 1 ,k 2 ,m)) 



(A2) 



where Sid(M) is the entanglement entropy of a one 
dimensional massive scalar theory with mass AI = 
^m 2 + 7 2 (2 - cos{ki) - cos(k 2 )). UsingS 2 . S 1D (M) cx 
-log(M 2 ), 



S a _ l °9( m2 + 7 2 ( 2 - cosih) - cos(fc 2 ))) (A3) 

fel,fe2 

Using Eulcr-Maclaurin formula, one finds the following 
expression for S*, correct to O(L ): 



S oc Ji + 1 2 



(A4) 



where 



h ~ L 2 log{m 2 ) +L 2 dt log(l + 7 2 /m 2 (l - cos(27rf))) 



and 



When rotational symmetry is broken, all powers of L 
(after the area-law term) are present, except in two di- 
mensions. In 2-dimensions, there is no constant term 
when the boundary of the region is smooth. So if rota- 
tional symmetry is broken, it is not as easy to recognize 
a topological phase in even dimensions higher than two 
without resorting to the ABC construction (Eqn. [15]) or 
its analogs. 

In one dimension, when symmetry is broken, it is con- 
venient to express the entropy in terms of x(s) and y(s), 
the parametric equation for the boundary. Because rota- 
tional symmetry is broken, there is no requirement that 
these terms appear symmetrically. However, because of 
translational symmetry, the entropy is a function only 
of the derivatives of these functions. The only require- 
ment is that the entropy density must be symmetric un- 
der s — > — s. Otherwise the entropy depends on whether 
s is measured clockwise or counterclockwise around the 
region, and this violates the symmetry between the inside 
and outside. (Clockwise is defined relative to a choice of 
the region's inside, as is familiar from using residues to 
evaluate integrals in the complex plane.) The expression 
J(^f) 3 rfs is symmetric and it scales as j?. One can 
check that the integral is nonzero for an ellipse (assume 
the ellipse has a small eccentricity so that the integral 
can be evaluated-one can then expand it to first order in 
the eccentricity). 

To see that there is no scale-independent term, notice 
that such a term would have to result from integrating 
an entropy density with units of 1/L. Such terms are of 
the form: 



F(s) = /( 



dx 
ds 



dy,.d 2 x dx 

Ts ){ d^ ) + 9{ Ts> 



dy, 
ds' 



ds 2 ' 



(Bl) 
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where s is the arc-length. 

These terms are total derivatives: let a define the 
angle of the tangent vector. Then ^| = cos a,-^- = 
sin a. The term equals [—/(cos a, sin a) sin a + 
<7(cos a, sin a) cos a] ^ . 

Now the total entropy, J F(s)ds can be rewritten as 
an integral with respect to a: 

J [—/(cos a, sin a) sin a + g(cos a, sin a) cos a]da (B2) 

By the symmetry s — > — s. / and g must be even functions 
of cosine and sine. Therefore the integral is equal to 
zero since the contributions from a and a + ir cancel one 
another. 

While there are no constant terms for smooth regions 
in two dimensions, shapes with corners do have scale- 
independent terms that can be attributed to the corners 
(this can happen even if the rotational symmetry is not 
broken, as shown for the quantum Hall state in Ref^i ) 

The first anomalous term in two dimensions scales as 
-p- ; one example of such a term is 



ds 



d 2 x 



(B3) 



This expression does not vanish identically, as can be 
seen by an example of a region which is nearly circular, 
i.e., described by the polar coordinates 



r = l + e(6), 



(B4) 



where eCl. Evaluating the integral to linear order in e 
one obtains 



ds 



ds 2 



ds(cos 3 0)(1 - 6e(0) - 3e'(0) tanf 



(B5) 

which is nonzero for e(0) oc cos 30. 

In higher dimensions, the entropy can have terms that 
depend on and higher derivatives. When rotational 
symmetry is broken, the Cartesian indices do not have to 
be contracted but the indices for the coordinates on the 
surface still do, because the entropy has to be indepen- 
dent of the coordinate system. The first allowed correc- 
tion to the entropy scales as — j — and a simple term in 
the entropy density that leads to such a correction 



Du a DuP / \du-r du gy 



(B6) 



The first factor gives the proper scaling, The other 
factors ensure that the expression is not a total derivative 
in dimensions above two. (The simplest term that scales 

V 2 x is a total derivative.) The 



as i — 

L ' Di 



' D 



j?9 



a/3 



indices are all contracted in pairs so the expression is 
coordinate- invariant. 

Once a term of order L D ~ 2 has appeared, one would 
expect that all terms of lower orders appear too, and 



that is what one finds. Since the L° term is especially 
important, we checked explicitly that in even dimensions 
greater than two, a scale-invariant term L° is allowed 
as long as there is no symmetry. As an example, take 
Eq. (|B6|) and multiply it by a power of the mean curva- 
ture, (g a ^ ' n a p) D 2 . This has the units L 1 ^ , therefore 
it gives a scale-invariant contribution when integrated 
(note that since D is even, there are an even number 
of powers of k as required). Both Eq. (|B6j) and this 
term with the extra factors of k give a nonzero entropy 
for a generic region. As an example, consider the D- 
dimcnsional surface of revolution obtained by rotating 
Eq. (|B4j) (which can be regarded as a curve in the xi, x^ 
plane in _D-dimcnsions) around the a;i-axis (explicitly, 
x x = r(0) cos 0, v / (x 2 ) 2 + ■■■ + (x D ) 2 = r{6) sin 0). Both 
integrals are nonzero for e(0) oc cos 30. 



2. Broken Parity Symmetry 

Breaking just parity can also lead to terms whose ex- 
ponents deviate from Eqnf4]even if rotational symmetry 
still exists, when the dimension is odd. However, these 
additional terms all vary as a negative power of L. In 
even dimensions, parity does not lead to any additional 
terms. 

The extra terms arise from an additional tensor, 7, that 
is allowed when parity is broken. This tensor is related 
to the antisymmetric tensor e: 



1 



7' 



V9 



c a±a>2 ...a>n-i 



(B7) 



where y/g is the square root of the determinant of the 
metric. (This coefficient is necessary for ensuring that 7 
transforms as a tensor). 

To understand these results, first note that the 7 ten- 
sor, like the n tensor, depends on how one chooses the 
normal to the surface. The sign of 7 depends on how 
the orientation of the surface is chosen, and this in turn 
depends not only on the orientation of space (which is de- 
termined by the parity-violating ground state) but also 
on the normal to the surface, see figure [4] Since 7 is odd 
under changing the sign of n, if a factor of 7 appears in 
the entropy-density, an odd number n K of factors of K 
must appear as well 



1( mod 2) 



(B8) 



Now the requirement that all the indices of the k's and 
its derivatives can be contracted with the upper indices 
of the g's and the factor of the 7 still implies that no, 
the number of covariant derivatives is even if D is odd, 
because then 7 has an even number of upper indices. So 
a term that includes a factor of 7 has units of L~ nK ~ nD 
which is an odd power of 4-. 

In an even number of dimensions, the entropy still goes 
down by two powers of L at a time because in this case, 
there must be an odd number of derivatives as well as an 
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Figure 4: Defining the orientation of a hypersurface 
from an orientation of space, illustrated in three dimen- 
sions. A pair of axes on the surface a, b is defined 
to be right-handed if the triad a, b, h is right-handed. 
Formally speaking, 7 is defined by contracting the D- 
dimensional epsilon tensor with the normal h and then 
transforming to curvilinear coordinates, 7° 1 ° 2 "" 1 '" 1 = 

n e ni2-.-i d Q^WT ■ ■ ■ du U D -i 9 9 - 9 



odd number of factors of ft to respect both the rotational 
symmetry and the Zi symmetry between the inside and 
the outside of the region. 

The first anomalous term in the entropy is however 
very small, and scales as or depending on whether 
the dimension of space is 1 or 3 modulo 4 respectively. 
All the terms between the area law term, L 15-1 and the 
constant term go in steps of L 2 . This is essentially be- 
cause the anomalous terms include factors of 7 which has 
D — 1 upper indices which all need to be contracted with 
something, forcing the term to have at least D factors of 
« or covariant derivatives, (one might think that y fac- 
tors of ft should be enough since each ft has two indices; 
however, since ft is symmetric and 7 is antisymmetric, 
contracting both indices of k with 7 gives 0.) 

For illustration, here are some examples of non-zero 
terms: 



h = J dAg aT 7 a ^D ai K a2a d T (tv ft) 2 

^5 = / dAg 1 "7 1 2 3 4 Z) ai k Q2(Ti Z) Q3 ft Q4Cr2 tr k 

for 3 and 5 dimensions, where tr k is the mean curvature 
times D — 1, .g Q ^K Q ^. These both scale as j^. Both 
expressions can be generalized to higher dimensions by 
introducing extra factors of D ak n ak+ia and adding 

factors of g" i<T i to contract all the cr's. 



Appendix C: Valid Constructions to Extract 
Topological entropy in D=3 

Let us define D[X(A, B,C)} where X denotes some 
property of a manifolds A,B,C to be 

D[X(A,B,C)} = X a +X b +Xc-X ab -X bc -Xca+X A bc 

For example, in this notation jtopo = D[S(A, B,C). 
We note that there is one obvious constraint for the 
Eqn. [15] to be useful which is that the number 
D[b Q (dA,dB,dC)} ^ 0. This is because 
D[bo(dA,dB,dC)} where the proportionality constant is 
the universal topological constant associated with the 
phase of matter. 

We claim that in D = 3, the Eqn. [T5]still holds as long 
as the following condition is satisfied: 

D[x(dA,dB,dC)} = (CI) 

where x denotes the Euler characteristic. This is be- 
cause this Eqn. IC1I guarantees that any dependence on 
local curvature cancels out on the right hand side of Eqn. 
[T5l Since \ = 2 — 2g, the above Eqn. may be re-expressed 
as D[g{dA, dB, dC)] = 1. Just to illustrate this point, 
consider the case when AnBnC 0. Generically, three 
regions in D = 3 would intersect along a line. Therefore, 
there are two possibilities: the region C wraps around 
the surface defined by An B or it doesn't. First consider 
the former possibility. In this case, as one may easily 
check that 

g(d(AC)) = g(dA) + g(dC)-l 
g(d(BC)) = g(dB)+g(dC)-l 
g{d(ABC)) = g{d(AB))+g(dC)-l 

This implies that D[g(dA,dB,dC)] = 1 and therefore 
this is a valid construction. An example is provided by 
the construction in Figj3^,. On the other hand, when C 
does not wrap around An B, then 

g(d(AC)) = g(dA)+g(dC) 
g(d(BCj) = g(dB)+g(dC) 
g(d(ABC)) = g(d(AB)) + g(dC) 

and therefore Q[g(dA, dB, dC)] = which implies that 
this is an invalid construction. 



Appendix D: Entanglement Entropy of layered Z2 
topological phases 

Here we briefly discuss the layered topologically or- 
dered states mentioned at the beginning of Sec. [TTJ These 
phases would lead to a correction ASU = —J2dL z , for 
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C 


B 





(a) 



(b) 



Figure 5: Extracting entanglement entropy when a three di- 
mensional topological ordered state coexists with a layered 
two dimensional topological order. Fig. (a): The ABC con- 
struction for this geometry yields 7 = L z "/2d- Fig.(b): The 
ABC construction for this geometry yields 7 = 73d . 



layering perpendicular to the z direction where 72D is 
the topological entanglement entropy associated with the 
theory living in each layer. For a generic geometry, the 
total topological entanglement entropy 7 may be written 
as 




Figure 6: Illustration for proving that the entanglement en- 
tropy is linear in the first Betti number, a) shows a genus 
k + 1 torus (k = 3), and b) and c) give two ways of dividing 
it up into A, B,C regions. In b) A is a fc-torus, B is a ball, 
and C is also a ball that has been stretched out. In c) A 
is a fc-torus, B a 1-torus, and C a ball filling up the hole of 
the 1-torus. Applying strong subadditivity to both configura- 
tions gives a recursion formula for S(k + 1) in terms of S(k). 
Strong subadditivity gives inequalities, but the left and right 
hand sides of the inequality in (b) and (c) are reverses of one 
another, leading to an exact relation. 



j{L z ) = 2j 2D L z +73D (Dl) 

Here L z is the dimension of regions A, B, C in the z 
direction. Fig |5] shows two different geometries for which 
the application of ABC formula (Eqn. IT5j) yields 7213 
and 73 d separately. 



Appendix E: Linear dependence of Stopo on Betti 
numbers in Three Dimensions 

Recapitulating the results from the previous section, 
in D = 3, the topological part of the entanglement en- 
tropy is proportional to bo and is independent of b\ since 
the genus dependence can be obtained by patching lo- 
cal Gaussian curvature. Implicit in these statements is 
an important assumption which is that entanglement en- 
tropy depends linearly on bo, b\. This form for the en- 
tropy also assumes that the entropy does not depend on 
knotting or linking of the toroidal surfaces. In this sec- 
tion we provide a proof of these two statements. Our 
only assumption is that the space in which region A is 
embedded is flat (i.e. has the topology of K 3 ). 

We begin by proving the linearity of the entropy for 
the simplest geometries shown in Figs|Bk and [7^,. Fig(SJt 
can be used to show that the entropy is linear in b\ . Con- 
sider first the decomposition shown in FigjGp where the 
boundary of region A is a k— torus (i.e. a torus with 
genus k), B and C are three-balls B 3 that join together 
so that the boundary of region B U C is a torus with 



genus one. We denote by S topo {k) the topological part 
of the entanglement entropy corresponding to a region 
whose boundary has genus k. Now apply the strong sub- 
additivity inequality, S AllB + S B uc > S A ubuc + S B - 
The local parts of these entropies satisfy an equality, 
Si oca i(AuB)+S local (BUC) = Si oca i{AuBUC)+S local (B) 
for the configuration shown, and any configuration where 
A and C do not meet. The reason is that each patch on 
the boundaries of the regions occurs an equal number of 
times on both sides of the equation. Hence the topological 
parts of the entropy also satisfy 

Stopo(k) + Stopo(l) > Sto P o{k + 1) + S tO po(0) (El) 

Similarly, we construct a different geometry as shown 
in FiglSJ; where the region A is the same as in FigJBp, 
region B has a boundary with genus one, and region C 
is topologically a ball B 3 . For this topology, adding C to 
AU B decreases the number of handles from k + 1 to k, 
so the strong subaddivity inequality yields the opposite 
conclusion: 



Sto P o(k + 1) + S topo {0) > S topo {k) + S topo (l) (E2) 

The above two equations imply that St opo (k + 1) — 
Sto P o{k) is independent of k and hence St opo is linear 
in k. This genus dependence can be traded for a non- 
topological contribution plus a 60 -dependence. 

To prove that the topological term is linear in the Betti 
number bo (= number of connected components of the 
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(a) 





Figure 7: Linear Dependence of the Betti number on 60. a) 
A spherical region with k + 1 hollow cavities in it (60 = k + 2 
including the outer surface). b,c)Two sets of regions for ap- 
plying strong subadditivity in order to prove that the entropy 
for k + 1-holes minus the entropy for k holes is a constant, 
b) The region is divided up into C, containing k — 1 of the 
cavities, a slab B, and a region A containing the remaining 
cavity, c) C again contains k — 1 of the holes, B is the rest of 
the region, and A is a solid sphere that fills one of the cavities. 



boundary) we repeat the above argument by replacing k- 
tori with a region with k spherical cavities cut in it (see 
FigE]). This yields the result that S top o = abo + (3- Since 
Stopo should vanish for 60 = 0, this suggests that St opo is 
strictly proportional to 60 . 

Now we will show that the entropy of any three di- 
mensional region, no matter how knotted, is given by 
the same formula in terms of the Betti numbers. First, 
we assume that 60 = 1. The theory of surfaces shows 
that every connected surface embedded in three space is 
topologically equivalent to one of the fc-holed tori. This 
torus may be built up by attaching handles repeatedly 
to increase the genus. However, there are many ways to 
attach a handle. For example, one could add a knotted 
handle, like in Fig(5^,. The argument in Fig(5J: does not 
apply to this handle because there is no way to define 
region C- for the argument above to work it should fill 
in the hole in the handle so that the genus decreases by 
1, but this is not possible since the handle is knotted 



Figure 8: Illustration of gluing and drilling, a) Gluing on a 
general handle. The torus region is Rk, and the shaded region 
is H, the handle. The handle meets the boundary of region R 
in a disk, which is shaded in black. Attaching the handle gives 
Rk+i = RkUH. b)Drilling a hole The initial region here, Rk, 
is a solid sphere, and H is a knot that is hollowed out from 
the inside of the ball to obtain Rk+i — Rk — H. The boundary 
of this region is topologically a genus one torus, although it is 
not possible to continuously deform it into one. Hence drilling 
increases the genus by one just as adding a handle does. 



through itself. 

To generalize the argument we therefore construct a 
proof of the upper bound that does not require filling 
the hole in. We will choose regions A,B,C that are all 
subregions of the initial region and the handle itself, so 
that the argument works if there is some linking. 

Ways of Building up a region First we will describe all 
the ways a region whose boundary is a connected genus 
k surface can be built up starting from a sphere Rq, so 
we can be sure that our argument applies to all of them. 
Given a region R k whose boundary is a fc-torus, one may 
either glue on a handle H as illustrated in FigJS^, or drill 
out a hole, as in Fig(5jD. A handle is a solid cylinder out- 
side the region R k whose two ends are on the boundary 
of Rk- Attaching H means defining Rk+i = Rk U H, 

At some stages it may be necessary to drill a hole out of 
the region instead. Drilling a hole is an inverted version 
of the previous process (see Figj8p): let H also have the 
topology of a cylinder, but H c Rk, with its circular 
faces on the surface of Rk- Then drilling the hole H 
is passing from Rk to Rk+i = Rk — H. For example, 
applying this process to a ball, as in Fig[8]3 gives a "knot 
complement," a manifold that is very interesting to knot 
theorists, since it encodes the structure of the knot. 

We will now prove that S topo {.Rk+i) ~ S t0 p {Rk) = 
St op o(torus) — Stopo(sphere) for the general case. We will 
focus on the case where Rk+i is constructed by adding 
a handle. (The other case is similar: if Rk+i is obtained 
by drilling, then the complements, R% +1 and are re- 
lated by handle- adding, so we can apply our arguments 
to these complements instead.) 

Change of entropy on adding a handle not affected by 
links passing through H We will first show that the en- 
tropy changes by Storus — Sbaii even if the handle H is 
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Figure 9: Figure for the proof that links passing through the 
handle do not affect the entropy. Part (a) shows the region Rk 
and the handle (in gray) that is added to it to get Rk+i- As 
in Fig[6l the goal is to relate the entropy of a k + 1-holed torus 
to a fc-holed one, where k = 3. The dotted circle indicates 
the portion of the figure that is enlarged in the later frames 
of this figure. The portion of Rk that passes through the 
hole is not shown in the subsequent frames for clarity. Parts 
(b) and (c) show the two decompositions that are used to 
prove the lower bound on Stopa(Rk+i) — S top o{Rk)- In (b), 
there are three regions: Region B (checkered) is Rk- The 
handle is flattened so that it is ribbon like, and then strip- 
like regions A and C are demarcated along its edges. These 
together form Au B U C := Rdoubie which has the topology 
of Rk with two handles attached. In (c) we decompose Rk+i 
into A, the center of the ribbon, B (checkered), the border 
of the ribbon closed up with parts of Rk to form a loop, and 
C, the rest of Rk- In this figure B U C = Rdoubie- Part 
(D) shows the decomposition used to prove the upper bound, 
which corresponds to the division in Fig[6b. 



linked with other portions of Rk- In this argument, we 
assume that the handle is not knotted with itself; the 
next argument shows that knotting does not affect the 
entropy. 

Consider the regions in FigJSJi. We will first prove this 



inequality: 

Stopo(Rk+l) — Stopo(Rk) > Storus ~ Shall (E3) 

where Storus and Shall are the topological entropies of an 
unknotted torus and a ball respectively. Begin by sliding 
the two ends of the handle (the gray region in FigEK) 
long the surface of Rk so that they are right next to each 
other, as shown in the figure. The proof of this inequality 
has two steps, using the decompositions illustrated in 
FigEp,c. 

First, as in Fig[9jD, let B = Rk, and let A and C be 
narrow strips along the left and right side of H. A U B 
and CUB are two regions that can both be deformed into 
the same topology (Rk+i); altogether Ai) B U C forms 
a region Rdoubie with two handles attached to Rk- The 
strong subadditivity implies 

2<5 topo (i?fc + i) > Stopo(Rk) + S t opo{Rdouble)- (E4) 

Next, as in FigJSJ:, let A be the interior of the strip 
H . Let B be the border of this strip, marked with a 
checkered pattern in the figure. This border is closed 
by adding some small parts of Rk so that it surrounds 
region A completely. Let C be the rest of Rk- In this 
decomposition, CUB = Rdoubie, while AUBUC = Rk+i, 
and AuB and B are a ball and a torus respectively. Then 

Stopo(ball) + Stopo{Rdouble)) > Storus + Stopo(Rk)- (E5) 

Adding the previous two equations and canceling Rdoubie 
gives Eq. (|E3|) . 

In this argument it has not been necessary to add on re- 
gions external to the handle, so the argument still works 
if parts of Rk are linked through it. 

Now the reverse inequality S t0 po(Rk+i) < S t0 po(Rk) + 
Storus — Shall is proved just as in the argument for the 
simple /c-tori using the division shown in FigEJd, which 
is no different than the original construction in FigJBb. 
Hence S t0 po(Rk+i) — S t0 po(Rk) = S tor us — Shall, and we 
conclude that S t0 po(Rk) = kS t0 rus ~ {k— l)S ba ii; that is, 
the entropy is linear in the genus of the surface. 

For the case where the surface has more than one 
boundary component we just start with a region whose 
boundary has many components, which are all spheres. 
The entropy of this region is proportional to bo by an 
argument similar to the one illustrated in Fig[7] (this ap- 
plies even when the spheres are nested in each other). 
It is known that any region may be built up from such 
a region by either attaching handles or drilling holes. It 
is even possible to choose the handles and holes so that 
each one starts and ends on the same component of the 
boundary of Rk ■ Then the calculation of the entropy as 
the handles are attached proceeds just as above. 

Change of entropy not affected by knotting of the han- 
dle If the handle is knotted, this argument does not im- 
mediately work, because the entropy change Storus — Shall 
could be different when the torus is knotted. However, 
a knotted handle may be transformed to an unknotted 
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Figure 10: Showing that the entropy of a knotted handle is 
the same as if it were an unknotted, (a) Shows the knotted 
handle H (the gray region) that is attached to the region Rk 
(the white region). Passing one strand through the other so 
the topology changes to figure (e) causes the handle to become 
unknotted: the loops in (e) can be untwisted, giving a simple 
handle. Any knot can become unknot-able if the right strands 
are passed through each other. The intermediate panels show 
that this process does not change the entropy. Panel (b) is a 
blow-up of Fig. (a). In panel (c) the topology is changed so 
that one strand passes through an eyelet in the other. This 
configuration is obtained by adding a handle along the dotted 
lines in (b), hence the entropy changes by Stoma — S Bp here- 
(The handle is unknotted) Panel (d) is obtained by removing 
a handle from (c)-the entropy therefore decreases back to its 
original value. Thus, the change from (a) to (e) does not 
affect the topological entropy. 



curve. This formula can be understood by introducing 
a new coordinate system where the x-axis is tangent to 
the surface at Pq. Then the y- value of the curve has 
a minimum at x — and the curvature k, intuitively, 
should describe how rapidly the curve curves away from 
the x-axis, i.e., y = H(x) = inx2 + ■ ■ ■ , which leads to 
Eq. IFT1 because H(x) is given by — h(P ) ■ (R(x) — R(x )) 
where x is used to parametrize the curve. (Although 
the derivative in Eq. (|F1[) is taken with respect to s 
rather than x, the answer is the same because x{P) rj 
s(P) — s(Pq) because the x-axis is nearly coincident with 
the curve in the vicinity of Pq.) 

In higher dimensions there are more directions avail- 
able for the surface to curve in, so the curvature is given 
by a tensor. Again one begins by rotating the coordi- 
nate system of space again so that it is easy to describe 
the curvature of the surface and then rewrites the answer 
in terms of curvilinear coordinates on the surface itself. 
Take the first D — 1 coordinates to span a hyperplane 
tangent to the surface, and let the D th coordinate point 
opposite to the normal vector; the surface is described 
by 



D-l 



H(xt 



1 



,XD-1, 



£2 



(F2) 



It is convenient not to have to keep redefining the co- 
ordinate system for each point one studies. One can 
define k instead relative to a fixed curvilinear coordi- 
nate system on the surface, described by the coordinates 



v},u 2 ,...u D 1 . Wc obtain 



one by repeatedly adding and then removing of unknotted 
handles as illustrated in FigflOl By the previous argu- 
ment, the entropy returns to its original value after this 
is done. Hence the knotted and unknotted handles have 
the same entropy. 



Appendix F: A lightening primer on differential 
geometry 

We have assumed that the entanglement entropy of a 
region can be written as a sum of contributions along 
the boundary. These contributions can depend on the 
shape of the boundary; for a smooth boundary, the shape 
is described, at lowest order, by the curvature tensor. 
The purpose of this section is to describe the curvature 
tensor and some of the operations that one can perform 
on this tensor to get the quantities that the entanglement 
entropy can depend upon. 

In two-dimensions, the boundary is a curve, and the 
curvature of such a curve at a point Pq is given by 



n(P ) 



d 2 R 

ds 2 



(Fl) 



where h is the normal vector at the point, R is a vector 
tracing out the curve and s is the arc-length along the 



d 2 H 



dx 1 dx° 
13 du a duP 



du 13 

d 2 R 

du a dv,P ' 



(F3) 



The first equation is the transformation rule for covariant 
tensors and the second uses Eq. (|F2j) and the chain rule. 

In three dimensions, the two "principle curvatures" k± 
and K2 are defined to be the eigenvalues of Kij relative to 
the special Cartesian coordinate system. They are usu- 
ally defined more geometrically: let a plane be drawn 
through the normal and allow it to sweep through all 
orientations. For each orientation, one obtains an inter- 
section curve and the largest and smallest curvatures of 
these curves are the principle curvatures. One can set 
up a special coordinate system aligned with the principle 
directions since they are orthogonal to one another. The 
two scalars that are most useful are the mean curvature 
and the Gaussian curvature, defined as + K2) and 
K\K2- These expressions can be written for any Cartesian 
coordinate system as ^tr n and det k. Generalizing these 
expressions to curvilinear coordinates cannot be done 
without introducing the "metric." These expressions do 
not apply to the curvilinear coordinates, though, since 
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'cur v/lt i tea r 



Mn Cartesian M T where M = |f . There- 
fore tr ^curvilinear = tlM Kc artesian M T , which does not 

simplify unless M is special. 

The technique for making coordinate-invariant expres- 
sions is to multiply tensors with equal numbers of upper 
and lower indices and to contract them all. The entropy, 
for example, cannot depend on the coordinate system 
because only the surface itself is important, not how it 
is described. Defining the mean curvature and Gaus- 
sian curvatures for a general coordinate system is a good 
illustration of how one makes coordinate-invariant ex- 
pressions. For any D — 1 dimensional surface there are 
always two tensors that one can define. The curvature is 
the more complicated of these, and the other is the met- 
ric which is the matrix of inner products of the tangent 
vectors along the coordinate directions: 



9a/3 



<9R 



(F4) 



which describes how stretched or sheared the u coordi- 
nates are. If one does not know the metric g, the compo- 
nents of K a /3 do not mean anything. They might be very 
big, but this could be because the mesh spacing in the 
it-coordinate system is gigantic. Instead, one can choose 
a coordinate system that is orthonormal according to the 
metric. More conveniently, one usually defines g a @ to be 
the inverse of the metric tensor and then contracts it with 
k in various ways, to get scalars For example, define 



.4 



D 
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(F5) 



These arc both independent of the coordinate system, 
and they can be calculated by using the orthonormal co- 
ordinate system with the axes along the principle direc- 
tions. A and B reduce to k\ + K2 and k\ + Hence 
H = iA, G = ^(A 2 — B). (Rationale for these calcu- 
lations: The tensor C" = g a ^ np-y maps vectors to other 
vectors of the same type, v a —> C"t> 7 (unlike the orig- 
inal k), so its eigenvalues arc well-defined. One cannot 
define the eigenvectors of a matrix that maps one vector 
space to another, because there is no way to compare 
the vectors to see whether they are proportional to one 
another.) 



The Gauss-Bonnet theorem states, interestingly, that 
the integral of the Gaussian curvature is a topological 
constant, equal to 2ir times the Euler characteristic of 
the surface; Chern generalized this result to manifolds of 
any even dimension: 



M 



d d 1 AAct k = -Lo d x(M) 



(F6) 



where d is odd, uid is the area of a sphere in (i-dimensional 
space and x(M) is the Euler characteristic of the man- 
ifold. (This result applies to d — 1 manifolds that need 
more than d dimensions to be represented without self- 
intersections-it only has to be rewritten in terms of the 
intrinsic curvature R a py§.) 

The last important tool from differential geometry is 
the covariant derivative; it is a generalization of gradi- 
ents from scalars to vectors. The most straightforward 

way to define a derivative of a vector does not have 
any meaning-it measures not just the variation of the 
vector, but also the variation of the coordinate system. 
For example, a vector field that does not change along a 
certain direction may have changing components because 
the coordinate axes may bend. To measure just the vari- 
ation of the vector, one can first transform to coordinates 
that are not rotating-namely, the Cartesian coordinates 
of the space, and then differentiate the vector. After this, 
one projects the change in the vector onto the surface, 
since the component of the derivative along the normal 
is determined by the fact that the surface is curving. In 
equations, one defines 



d-i 



W = J2 




7737 • (F7) 



Although this definition does not make it clear, the co- 
variant derivative can be calculated without using a co- 
ordinate system for the embedding space; it can be ex- 
pressed entirely in terms of the metric g a /3, as is familiar 
from general relativity and differential geometry books. 
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